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1 Introduction 

The program of Stochastic Deformation was born in 1984-5 as an attempt to un- 
derstand the paradoxical probabihstic structures involved in quantum mechanics 
[1] . In the course of this work, it became clear that no mathematically consis- 
tent and physically relevant approach was (and perhaps ever will be) available. 
But also, on the positive side, that a mathematical reinterpretation of Feyn- 
man's Path Integral approach to the problem was indeed accessible, providing a 
general method to produce quantum-like probability measures with qualitative 
properties quite distinct from what we are used to in statistical physics. 

It is our intention to describe here the main features of this method. No 
familiarity with Feynman path integral approach or even quantum mechanics 
itself are required. We shall summarize their basic elements allowing the reader 
to understand why our construction is a rigorous version of his approach. 

Quantum Theory provides for sure a kind of deformation of Classical Me- 
chanics. But not a probabilistic one, unfortunately. One way to see our con- 
struction is, precisely, as such a quantum-like stochastic deformation. In this 
respect, it should also be of interest in stochastic approaches to Geometric Me- 
chanics [2]. 

It may be worth recalling that the project to make sense of Feynman Path 
Integral method (and not only of a few of his formulas) is still a widely open 
problem, of a potentially devastating generality for probabilists. Who can doubt, 
indeed, after all these years and so many applications far beyond what Feynman 
could imagine in the fifties that a mathematically consistent version of this 
method should exist ? 

This overview of Stochastic Deformation will be organized as follows: 
Section 2 will summarize the original ideas of Feynman's reinterpretation 
of elementary quantum mechanics. For a probabilist, they look like a (very) 
informal version of Stochastic Analysis, with a twist regarding boundary con- 
ditions. This twist will prove to be indeed fundamental for their mathematical 
interpretation. 

Section 3 will provide a probabilistic counterpart of Feynman's approach, 
i.e. a kind of stochastic boundary value problem whose conditions of existence 
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and uniqueness of solutions will be specified. 

In Section 4 the random dynamics of the relevant class of processes, together 
with its associated symmetries, will be described, in Lagrangian and Hamilto- 

nian form. 

The last Section will be devoted to some computational and geometric as- 
pects of our Stochastic Deformation. 

It is a pleasure to thank the organizers of the stimulating semester program 
"Stochastic Analysis and Applications" at EPFL, 9/1/2012-30/6/2012, namely 
Robert Dalang, Marco Dozzi, Franco Flandoli and Francesco Russo. The active 
participation of the audience helped to improve the presentation of the material 
and is also warmly aknowledged. 



2 Overview of Feynman Path Integral approach 
to quantum theory 

We shall consider the example of a system of a single unit mass (charged) particle 
under the effect of a (bounded below) scalar potential V{q) and a smooth vector 
potential a(q). According to elementary quantum mechanics, the Hamiltonian 
observable of such a system is a densely defined self- adjoint operator (in L^(R^) 
for simplicity) 

H=^{P-a{Q)r + V{Q) (2.1) 

where the position and momentum observables are defined on appropriate dense 
domains as 

V'(g) ^ q^{q) 

(2.3) 

^(9) H- -ihVipiq) ^ ' 

The initial state i\) of the system evolves in an unitary way, Ut"^ = V't) = 
e~t*^, so that tpt solves 

ih^jjt = mt (2.4) 

at 

Any observable, in fact, is a densely defined self-adjoint operator O. The 
contact with (the language of) probability theory is made by the definition of 
the "mean value of O in the state tjjt" at time t: 

(O)^, = (AlOiJt) = (V-IOtt/-) = {Ot)^, Ot = U+OUt (2.5) 

where (•!•) denotes the scalar product, antilinear on the left (with associated 
norm || • ||), and -|- is the adjoint. 
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Eq (|2.5p is a dual expression of the dynamics of an observable. For H as in 
(|2.ip . we obtain 

\4i{Pt~a)^^{{Pt~at)Amtat~yV{Qt))^ 

But of what kind of random variables, really, those expressions are the mean 
values? We do not know; in fact no probability space has ever been defined in 
the first place. The only hint at probability theory, in quantum mechanics is Von 
Neumann's Axiom of "Quantum Static" according to which if one performs a 
measurement of an observable O for a system in state ip, the absolute probability 
to find a result ^ A G R is 

{E^iX))^ = WE^mr (2.7) 

where E'~'[\) denotes the spectral family of orthogonal projections of O. For 
instance, when O = Q, and any A in the interval [a, b], 

\\EQ{[a,bm'=f imi^dq (2.8) 

J[a,b] 



called Born interpretation of the wave function -0. In this sense Eqs (|2.5H2.6p are 
understandable, if not justified probabilistically. Even for the simplest quantum 
systems (Eq (|2.ip with a = V = O, i.e the "free case") there is no underlying 
concept of space-time trajectory. The justification of this prohibition lies in the 
uncertainty relation: 

For ip in the domains of Q and P, 

{QjPk - PkQj)^ = ihSjk, 1 ^ fc ^ 3 (2.9) 

is interpreted as the impossibility to localize experimentally the position and the 
momentum simultaneously, i.e to define a trajectory as in classical Hamiltonian 
mechanics. 

Feynman transformed qualitatively these shaky relations between Quantum 
Physics and Probability Theory [31. 

Let us consider a classical (conservative) system of Lagrangian L, in dynam- 
ical evolution on the time interval I — [s,u\. If its configuration at time t G I 
is denoted by uj{t), the Action functional of this system is defined by 

Sl[lo{-);u-s]= [ L{uj{t),uj{t))dt (2.10) 



For our system whose quantum Hamiltonian is (|2.ip . 

L{uj,uj) = ^\uj\^ ~V{oj) + a{uj)-6j. (2.11) 

Feynman starts from two states ip, ip, say at the final time u of I and rewrites 
(</'«|'0m) G C in terms of the unitary evolution kernel applied to ipg- 

tps{x){e'i^''"'^"){x, z)ip^{z) dx dz (2.12) 
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where denotes the complex conjugate. 

The key point is that he reinterprets this as the foUowing "Path Integral" : 

V's(x)e^'^^["(-)'"""l2?w^„(z)dxdz (2.13) 

In this symbohc expression, r2§'" means the set of continuous paths {w G 
C([s, u],IR^)|aj(s) — x,uj{u) — z} and Vuj plays the role of a "flat measure" 
ris^tsgu dui{t) on this path space. 

Consider any time s < t < u. Clearly can be regarded as initial condition 

of 

= H^, (2.14) 
and ^„ as final boundary condition of the adjoint problem 

- ih^ - Hip (2.15) 

Feynman calls "transition amplitude" the (complex) expression (|2.12p or 
(|2.13p . He will use it as a kind of expectation with weight exp j^Sl, and denotes 
it by ((ysjllT/i) or {1)sl- should be stressed that, in Feynman's view, (|2.13p 
should simply be regarded as a short notation for a sum over time-discretized 
paths uj{tj) = Xj, tj = j ''"j^^'' , 1 ^ j ^ iV G N. Along the same line, for "any" 
functional F[a;(-)] of such paths, the "expectation" of F is defined by 

{F)s,= I I /^s(a;) e*^-['^(-)-"-^lF[c^(.)]^Jz)Pwdxdz. (2.16) 



This definition is the starting point of Feynman's "Functional Calculus" , whose 
main result is an Integration by parts formula [3]: 

{5F[Lo]{5uj))s,^~'j^{F5SL[io]{5u:))s, (2.17) 

where 5G[ijj\{5uj) denotes the Gateaux derivative of a functional G on cj in 
the direction 5uj. Shaky as it is, mathematically, this formula deserves some 
interest since it is the ancestor of all the integration by parts formulas designed 
in Stochastic Analysis during the last 25 years ! 

For F = 1, L as in (12.111) . Feynman finds the path integral counterpart of 

{i;j)s^ ^{ujA rot aiuj) - \7V{u;))sr. (2-18) 

where, regarding the meaning of the l.h.s, he observed that "in the few examples 
with which we had experience, the substitution w = (^^^yi {^{t + ^t) ~ 2aj(i) + 
uj{t — At)) has been adequate". 

The most revealing application of the relation (|2.17p . however, is the one for 
F[uj] = uj and some specific direction 6uj. Indeed, (|2.17p reduces, then, to the 



4 



(discretized) expression: 



(2.19) 

This maybe the most fundamental discovery of Feynman approach, j3] although 
it is almost never cited, for reasons to be understood afterwards. First, let us 
observe that (|2.19p is a kinematical claim regarding the nature of trajectories. 
No dynamics is involved here. By construction, w(-) G ^^1'", but what (|2.19p 
suggests (since Feynman was aware that the desirable limAt4,o is, to say the least, 
problematic) is that w(i-) should be different from cj(<+), for any t £]s, u[. In 
other words, the a;(-) are indeed continuous quantum trajectories, but there are 
not differentiable. 

Now consider a = in the Lagrangian (|2.1ip . Then, by definition of the clas- 
sical momentum, p — ^ — uj ii the configuration q — lu. So Feynman regards 
(|2.19|) as the "space-time" version of the uncertainty relation p.9|) . providing 
a deeply different interpretation from the regular one of Quantum Theory in 
Hilbert space: no reference to any limitation on experimental measurement is 
involved here. 

Notice that, in this perspective and for Lagrangians as before, it would make 
perfect sense to interpret Quantum Mechanics as a Stochastic Deformation of 
Classical Mechanics for smooth paths. Clearly, h is the deformation parameter. 
However, even S'i[w(-)] would become singular along such "quantum paths" 
since lj does not make sense. 

Let us come back, for instance, on our question, after Eq (|2.6p : of what 
kind of random variables are those equations mean values? Take the above 
elementary momentum, for instance. By (|2.6p (with a = 0) it should satisfy, in 
particular, after quantization, 

= (Pt)^- (2.20) 
According to Feynman, however, this corresponds to 

, , / uj(t + At) -uj(t)\ 

where aj(-) denotes some, yet unspecified, random process. 

The right hand size can be computed as a difference of 2 functionals like 
(|2.16p . using Eqs (|2.14I - [2T5| . On the other hand, of course, according to regular 
quantum theory, and using ()2.3|) . we know that 

(Ft)v, - {P)i„ = I ^t{~^hVi>t)dq = /^t^t (-ih^\ dq. (2.22) 



In this way Feynman reinterprets each quantum observable as a specific func- 
tion of his "random process" uj{-), here —ihV logtpt for P. Unfortunately, R.H. 
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Cameron proved in 1960 ([!]) that e^^^Vuj does not make any sense as a count- 
ably additive measures on the path space underlying (|2.13l) . This means that 
there is no such thing as the "expectation" {)s^ of (|2.16l) . and no limAt4,o in 
expressions hke (PTTO)) or ^T^. 

A complex measure similar to the Wiener one but with a purely complex 
variance, like the one used informally by Feynman would be in particular of 
infinite total variation and therefore not appropriated for his quantum purpose. 

The challenge we face is therefore to preserve the essential of Feynman's 
approach but with well defined probability measures on path spaces. In order 
to tackle this, we will first summarize some of the key qualitative aspects to 
preserve: 

1. Assume that there is a (filtered) probability space (51,(7, P) where such a 
stochastic process, say X{t), t G I, is well defined. 

Then the uncertainty relation (|2.19p should mean 



where Et should be a conditional expectation since we know that the 
momentum has to become a function of the process X{t). This means, in 
particular, that X{t) should be Markovian. Still, the two limit random 
variables of ()2.23|) should be distinct, and will be denoted respectively by 
DfX and DtX. Following R.H. Cameron, we cannot hope to produce 
complex measures satisfying Eq (|2.19p . We would content ourselves with 
real measures. 

2. We need to preserve a probabilistic counterpart of the time symmetry in 
quantum physics (at least for systems whose Lagrangian is time indepen- 
dent) involved in Born interpretation (|2.8p : Indeed, the relation between 
ip and tp can be regarded as a time reversal (cf. (j2.14p - (|2.15|) ') therefore 
their product, in Born interpretation, is unchanged under this symmetry. 

3. Although we shall treat here exclusively the class of elementary Lagrangian 
systems (|2.1ip considered by Feynman, our construction should rely on 
general principles compatible with any physical system. 

Before describing the program of Stochastic Deformation, a few words about 
the internal evolution of Stochastic Analysis itself. 

After the pioneering works of Wiener and Ito, this field made great progress 
since 1980 but there is one where it did not; with the notable exception of some 
approaches of Stochastic Control (we shall come back on this), the field suffered 
from a chronic lack of dynamical content (in the classical sense of dynamical 
systems theory). In Ito's original perspective Stochastic differential equations 
are Stochastic deformation of ordinary differential equation, still the comparison 




(2.23) 
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of the history of the two fields is reveahng. For (second order) ODE a single, 
very hard, dynamical problem became the motor of all scientific progress: the N- 
body problem. Nothing like it was ever considered in Stochastic Analysis. This 
explains why very basic notions of ODE's theory, like the one of integrability, 
for instance, are lacking in Stochastic Analysis. 

Here is a (dynamical) "paradox" mentioned by Krzysztof Burdzy. Let : 
[0,T] — M such that sup^gjQ |^(i)| < oo. If Wt denotes the Brownian motion, 
it is known that the probability 

P{cf,{t) -e<Wt< m +e,Vte [0, T]} ^ c(e) exp -1 j\^{T)fdT = F[c^{-)] 

The functional F[<j)\ is maximized by (/'(t) = Vi G [0,r] i.e, in particular, the 
solution of the second order ODE 4){t) = with cj){0) = 4>{0) — 0. But what 
does it mean, really, in relation with any observation and known properties 
of Brownian paths ? Notice that, according to Feynman, this is the free case 
V = a = in (|2.1ip and the dynamical equation of this system should indeed 
be (cf (|2.18p ') {uj)sl = Oj whatever meaning can be given to {■)sl- 



3 Probabilistic counterpart of Feynman's 
approach 

To make sense of (|2.23p . in the form 



E 



Xi(t) lim Et 



X{t) - X{t - At) 



At 



lim Et 



k 

X{t + At) -X{t) 
At 



x,{t) 



= hS. 



•jk 



(3.1) 



we need two filtrations, to take into account not only the usual past information 
on a time interval / D [s, m], i.e an increasing one T^t, i G /, but also a decreasing 
filtration Tt taking into account the future. The underlying filtered probability 
space should, therefore, be of the form (17, cr, {Vt}, {J't}^ P) with t e [s, u]. 

Feynman's time discretized interpretation of the left hand side of his dy- 
namical Eq ()2.18p suggests to limit ourselves to processes Xt such that, for any 
bounded measurable / and any s ^ si < f < ii ^ u. 



E[f{Xt)\Vs,yJTt,]=E[J{Xt)\X{si),X{ti)] 



(3.2) 



This property is what we call now Local Markov (or two-sided Markov) . But it 
was introduced in 1932 by Sergei Bernstein, who named it "reciprocal" His 
motivation was a remark of E. Schrodinger, a year earlier, which seems to be at 
the origin of all the notion of stochastic reversibility known today to probabilist 
[B]. But we will refrain to insist here on the (quite tortuous) story of this notion. 

To keep track, as suggested by Feynman, of the past and future informations 
about the system, the traditional Markovian transition probability should be 
replaced by a more symmetric measure Q, named after Bernstein: 
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A I— >■ Q{s,x,t, A,v, z), \fx,z 6 R, s < i < It in /, measurable in x,y with 
Ae B{R), the Borehan of K (for simphcity). 

For Q there is a 3 points analogue of Chapman-Kolmogorov property, such 
that, for X{u) = z fixed, Q becomes a forward Markov transition and for 
X{s) = X fixed Q reduces to a backward Markov property (let us recall, with 
A.D. Wentzell 17^, that these are just two of the 64 ways to express Markov 
property !). 

Of course, without fixing the starting or ending point, X(-) will only be a 
Bernstein process, satisfying p.2p and not a Markovian one. 

Let us denote by M the joint probability measure on B x B for the pair of 
initial and final random variables. 

Then B. Jamison (1974 [8]) proved the following (cf. [1] or [17], for a version 
involving Vt and J-t and appropriate to the relation with Feynman's approach). 

Theorem. For a given Bernstein transition Q and a given joint measure M , 

a) There is a unique probability measure Pm such that under Pm, X{t), t G 
[s,u\, satisfies Bernstein property 



b) Pm{X{s) e As,X{u) e Au) = M{As x Au) Borelians in B, the Borel tribe 

c) Pm{X{s) e A,,X{ti) e Ai,...,X{tn) e A„,X(m) e A^) 
dM{x,z) / Q{s,x,ti,dxi,u,z) 

AsXAu J Ai J A2 



Q{tn~l,qn-l,tn, dqn, U, z) 



A„ 



for s ^ ti ^ t2 ^ ■ ■ ■ ^ tn < u and Ai E B, i = 1, . . . ,n. 

The final random variable has been fixed here so that, as said before, Q has 
the properties of a forward Markovian transition but c) would hold as well with 
a fixed initial X{s) = x. In other words, the construction is perfectly symmetric 
with respect to the past and future informations, as it should be. 

Jamison also proved that only one joint probability measure M = Mm {Mm 
for Markov) turn X{t) into a Markovian and not only a Bernstein process. 
Using the same notation as in §1 but, this time, for the (strongly continuous 
contraction) semigroup generated by the lower bounded operator H, Mm is of 
the form 

Mm{As xAu)= I iilix) (e-*("-^)^) {x,z)T^u{z)dxdz (3.3) 

where 77* (a;) and ri^iz) are two positive (not necessarily bounded) functions to 
be determined later. This expression should be compared with Eq (|2.12p . Now, 
with 

(5(s, x, dg, u, z) = h{x^ u — s, z)~^h{x, t — s, q)h{q, u — t, z)dq, s < t < u 
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where the handy notation (e~'^("~''^^)(x, z) = h{s,x,u,z) has been used, the 
substitution of p.3l) in the above finite dimensional distributions c) provides, 
after simphfications, the finite dimensional distributions: 

Pn{dxi,ti,dX2,t2, . . . , dXn,tn), S < ti < t2 < ■ ■ ■ < tn < U 

ri*{x)h{s, dx, ti, dxi) . . .h{tn, dxn,u, dz)rju{z) dx dz (3-4) 



Now define the following densities of the forward and backward transition prob- 
abilities 

P(s, x, dy) = h{s, X, t, y)^^^^^dy s ^ t (3.5) 

Vs[x) 

where 

rjsix) ^ / h{s,x,t,y)r]tiy)dy, (3.6) 



and 

P*{s,dy,t,x) = ^h*{t,x,s,y)dy s^t (3.7) 

where 

% (x) = J V*s {y)h* (t, x,s,y)dy^ J tjI iy)h{s, y, t, x) dy (3.8) 

and the classical relation h* (t, x, s, y) = h{s, y, t, x) between integrals kernels of 
two adjoint parabolic equations (for H not necessarily symmetric) has been used. 
Then it is easy to check that ()3.4p coincides with the final dimensional distribu- 
tions of a forward Markovian process of initial probability density rjl{x)rjs{x) dx 
and transition probability density p.Sp or, equivalently, of a backward Marko- 
vian with final probability density r]l^{z)r]u{z) dz and backward transition den- 
sity of the form p.7[) . 

As a matter of fact, Vi € [s,u\, it is true that 

PiXit) e A) - / vtVtix)dx (3.9) 

J A 

where, as shown by (|3.6[) and (|3.8p . rj^ and rjt are two positive solutions of the 
two adjoint parabolic PDE, s ^ t ^ u 

V [s,x) = ri^{x) 
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and 



r]{u,x) = riuix). 



(3.11) 



Here, a comment is needed since we said that quantum Hamiltonians, like 
any observable, are self-adjoint and the first PDE of (I3.10|) involves the adjoint 
ff+ of H. Our example (|2.ip illustrates this point. When written explicitly, 
Schrodinger equation (|2.4p means 



— AV + yV.aV■ 



(3.12) 



The "Euclidean version" of this corresponds to transform t into it and a into 
—iA so that the right hand side operator becomes, indeed, a non-symmetric 
operator 

H ^ -—A + hA\/ + -\/.A- -\A\^ + V. (3.13) 

Notice how close we are, in Eq p.9p . from Born interpretation of the wave func- 
tion ip in (|2.8p . Informally, we have done t i-^ it and the above transformation of 
the vector potential a. The product structure of the probability density in (|3.9p 
is fundamental; it expresses manifestly a kind of invariance under time reversal, 
more general than the one traditionally known by probabilists as "reversibility" . 
For instance, if Pt{dx) denotes this probability at time i, it follows immediately 
from (|3.5p and (|3.7p that the following "detailed balance" condition holds 



Ps{dx)P{s, X, u, z) = P* {s, dx, u, z)Pu{dz), s ^ u 



(3.14) 



generalizing Kolmogorov's notion of reversibility [9j to non-stationary situations. 
In this paper, in fact, Kolmogorov refers to E. Schrodinger (1931,32) [6] who is at 
the origin of our program of Stochastic Deformation. This has been regrettably 
forgotten afterwards. 

Given (|3.5p and (|3.7p . a simple calculation provides the forward and back- 
ward drifts of the underlying diffusion process X{-). Preserving the notations 
DtX and D^X introduced after jMSl), 



DtX = hm Et 

AtlQ 



D;X = lim Et 

AtiO 



X{t + At) -X{t) 



At 

X{t) -X{t- At) 
At 



hV\ogTit{X)-A{X) 
---hV\ogf^;{X)-A{X) 



In particular, since Pt{dx) = rj^rjiix) dx = p{x, t) dx 

D;X = DtX - Wlogp. 



(3.15) 



(3.15*) 



(3.16) 



After substitution of (|3.16p in the left hand side of p.ip and an integration by 
part we obtain an elementary proof of this uncertainty relation (j3.ip . justifying 
in this way the presence of two filtrations. 



10 



There is a new qualitative aspect in our probabilistic counterpart of Feyn- 
man's approach. Not surprisingly it comes from its boundary conditions. 

The above construction of the processes, for a bounded below "Hamiltonian" 
strongly suggests that natural boundary conditions should be two probability 
(densities) at the initial and final time: 



instead of the boundary conditions of the two adjoint equations p.lOp and 
(|3.1ip . But let us write the marginals of the joint probability Mm (|3.3p : 



Ivtix) J h{s,x,u, z)r]uiz) dz = P,(a;) 
[Vuiz) J r]l{x)h{s,x,u,z)dx ^ Puiz) 

in terms of the single integral kernel h. If Pg and P„ are arbitrarily given, Eq 
(|3.18p is a non linear integral system for {ri*,r]u) the two boundary conditions 
of the underlying adjoint PDE ((3T0l) and (|^TT|) . s < t < u. 
Beurling has proved in 1960 the following general result: 

Theorem ([IH])- Let the above integral kernel h{s, x, u, z) — (^^'ki'^^^)^^ (x, z) 
be continuous, positive and defined on any locally compact space. Then the 
system (|3.18p has a unique pair {ri*,rju) of positive, not necessarily integrable 
solutions, for any strictly positive probability densities Ps{x), Pu{z). 

The proof of Beurling uses an entropic argument. This approach is quite 
natural in many respects (cf. |llj for instance) when handling this class of 
processes. Here, however, our present motivation comes from Mechanics and we 
shall not elaborate this Statistical Mechanics connection. 

The above construction provides the solution of a stochastic boundary value 
problem, quite distinct from the Cauchy kind of problems originally inspired 
by Kolmogorov. Remarkably enough the processes solving such boundary value 
problems are necessarily invariant under time reversal (in a sense to be specified 
soon) although, as we will see, some of their partial characterizations reintroduce 
an "arrow of time" . 

Examples. 

1. Consider a Brownian Wt, t G R+, on the real line, with diffusion coeffi- 
cient h and initial probability density % > 0. This is the case A = V = 
in (j3.13p . i.e H — — -^A. The traditional interpretation is that, given 
li{dx) = x{x) dx, P^{Wt G dx) = ri*{x,t), where r/* solves 



Now, if we wish to look at Wt as a Bernstein reciprocal process X{t), we 
should start from the same Hamiltonian H , a bounded time interval, say 
I = [0,r] and the following boundary probability densities in (|3.18p ; 



Ps{dx) = Ps{x) dx and Pu{dz) = Pu{z) dz 



(3.17) 




(3.19) 
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Po{x) = and Pt[z) = r]^{x,T), where ry* is the (positive) solution of 
Eq (|3.19p . Of course, the kernel h of p.lSp is the Gaussian one: 

h{0,x,T,z) = (2^nT)-i/2cxp-|:t^ (3.20) 

and the solution of Eq p.lSp (the "Schrddinger system") on [0,T] is triv- 
ial: 

{7^*ix)=xix),VT{z) = l} (3.21) 

With those boundary condition, the solutions of the two (heat) Eqs. p.lOp 
are clearly, \ft € [0,T] 

ri*{x,t) ^r]^{x,t), r]{x,t) = l. 

According to p.isp the forward and backward drifts reduce therefore to 

DtX ^ 0, D;X = -Wlogr/*(X,t) 

so that, denoting by dX and d^_X , respectively, the ltd differentials under 
Et in ((3lB . X{t) solves both SDE 

dX{t) = h^/^dWt, d^X{t) = -hV\ogr^*^{X{t),t)dt + h^l'^d^W; (3.22) 

where denotes a Wiener process with respect to the filtration Tt, t ^ 
%T\. 

With such choice of boundary probability densities {Pq,Pt) it is hard, of 
course, to see any dynamical time symmetry. But let us switch them both, 
for the same kernel h as before. Then {?7o(a;) = l,fjx{x) — xi^)} '^^^^ 
solutions of the system (j3.18l) . Indeed, f}*{x,t) = 1 and 'q{x,t) = ri'^{x,T — 
t) solve the pair of heat equations p.lOp . p. lip . The new associated 
process X{t), t £ [0,T] is such that DtX = Wlogf)(X,t) and D;X = 1 
and it is as well defined as the above diffusion X{-). 

Notice that it follows easily from the above argument ( or directly from the 
definitions p.lSp and p.l5*P ) that 

DtX(t) = -D*X{u + s -t), ss^t^u. (3.23) 

This rule deforms the classical time reversal of derivatives into a more 
subtle one, involving necessarily two filtrations. 

The full time symmetry of a Bernstein process ( or equivalently of its prob- 
ability measure) appears more clearly when considering processes not of 
independent increments: 

2. For the same H as in 1. pick the informal limiting case Pg = and 

Pu — Sz, corresponding to the solution r]*{-) = Sx, rjui') — of Eq p. 181) . 

5*0 that rj*{q,t) = h{s,x,q,t) and r]{q,t) = h{q,t,u,z) with h as before. 
Eqs p.lSp provide the two drifts 

z-X(t) , Xit)-x 
DtX = D*tX = —-^ (3.24) 
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X{t) is called the Brownian Bridge between (s,x) and (u, z). Defining 
X{t) = X{u + s — t) , s ^ t ^ u, X{t) is another bridge traveling backward 
from Pu ^ 5z to P^^ 5^. 

3. The above construction is, in fact, independent of the form of the "Hamil- 
tonian" H. For instance, if one considers the (non-symmetric) 

Hri[k) = U{k)-q ~ cVrj - ^Ary- 

/ {ri{k + y)-r^{k)-yVr,{k)l{\y\^^yy{dy) (3.25) 

for U : M. continuous, bounded below, c and fc G M^, i^{dy) a Levy 

measure on ]R^\{0}, the resulting processes are well defined. They form 
an interesting class of time reversible Levy processes fUS^ . JJ^)- 

The above examples suggest the following general notion of the time re- 
versibility involved here: 

Any Bernstein process X{t), t e [s,u\ constructed as before, given a non 
necessarily symmetric Hamiltonian H with integral kernel h as in Beurling's 
Theorem, and any given pair of strictly positive probability densities Pg (x) and 
Pu{z) is invariant under time reversal in the sense that X{t) = X{u -\- s — t), 
t G [s, u] is also a well defined process of the same class, evolving backward from 
Pu to Ps. 

We shall conclude this section by a comparison between Bernstein measures 
and the usual ("Euclidean") approach in Mathematical Physics. Consider the 
definition p.3p of the Markovian joint probability measure Mm- For H self- 
adjoint and r/*,r]u real- valued and bounded it can also be regarded as the L^ 
scalar product 

and expressed in terms of Wiener measure /i^ : 

J 77:(c.(s))e-*^"^("(^))'^^774c.(zi))dAi^(c.). (3.26) 

This version of Feynman-Kac formula has been known and used since the sixties 
[57] ■ A key difference with our construction is that, to produce Bernstein mea- 
sures, 77* and rju have first to be found as (positive) solutions of the system (I3.18p . 
given initial and final probability densities Ps{dx) and Pu{dz). Only then the 
reversibility of Bernstein measures and therefore their dynamical meeting will 
show up. Indeed, in Eq (j3.26p the Wiener measure does not carry any specific 
dynamical meaning, in contrast with Bernstein measure, involving two drifts. 
This is already seen clearly in the Feynman-like Formula (3.1) and, of course, 
in the equations of motion that we are going to obtain in the next Section. 
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4 Stochastic Dynamics and Symmetries 



There are two approaches to classical dynamics (or the classical calculus of vari- 
ations) the Lagrangian and the Hamiltonian one. Feynman's method suggests 
to start from the former one. The classical system with Lagrangian (|2.1ip will 
be our guide. We already know that, in relation with the Vt filtration, the clas- 
sical Co should become DtX (apart from its imaginary unit factor) under our 
stochastic deformation. So, from the transformations used in (I3.12p . (13.131) the 
Lagrangian should be proportional to 

C{X, DtX) - i| A^P + V{X) + A-DtX + ^V -A. (4.1) 



The last term of (|4.ip requires some explanation , useful also for later purposes. 

In our time-symmetric context, 3 definitions stochastic of integrals are avail- 
able [14] for the stochastic deformation of a classical expression of the form 

A{u)duj{T) (4.2) 

Jt 

in the Action functional. 

With respect to the increasing filtration Vt, for A T^T-adapted, with 

E j^\A\^{X{T))dT < oo, 

/ A{X)dX{T)= l.i.p ^A(X(r,_i))(X(r,)-X(r,_i)) 

J t max \ Tj—Tj-i |— J-O - -^ 

where l.i.p means limit in probability. Then, using (|3.15p , 

AdX{T)^Ej ADrXdT. (4.3) 

With respect to the decreasing filtration 7> and introducing the notation d^X{T) 
for the backward differential involved in (|3.15*|) . 

/ A{X)d,X{r)^ l.i.p ^A(X(r,))(X(r,)-X(r,_i)) 

•It max |rj — Tj_i |— >0 



and therefore, under the J> counterpart of regularity conditions, 



/•U /'U 

E I Ad^X{T) ^E AD*XdT. 
Jt Jt 



(4.4) 



The third one is due to Stratonovich: (when, in addition, E J" IV A\^ {X {r jjdT < 
oo) 



/ AodX{T)= l.i.p ^-[A(X(r,_i)) 

i max |rj— Tj-i |— J-O - "i^ ^ 



A{X{TjmX{r,)^X{T,^,)) 
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and ^ ^ 

E AodX{T) ^ E A-^{DrX + D;x)dT. (4.5) 

It follows from the usual (forward) Ito calculus (TS] that 

AodX{T) ^ A-dX{T) + ■ Adr (4.6) 

Together with (|4.2p this justifies the following definition of the Action functional 
relevant for the dynamics of our system with Lagrangian (|4.ip : 

J[X] = E,,t ^ CiX{T), DrX{T)) dr 

= E^^t l\DrX{T)\^ + V{X{t)) dT + J\o rfX(r)| (4.7) 

where E^^t denotes the conditional expectation given X{t) — x, t < u. 

What about the geometrical meaning of J ? Let us define, for any C, the 
(forward) Momentum as classically by 

assuming that Eq (|4.8p is solvable in DtX ~ (p{P,X). This is the case for the 
system (|4.ip and 

DtX = P-A. (4.9) 

Then, if Qpc is the section, on the (X^t) submanifold, of the Poincare-Cartan 
1-form, the Action can be written as 



J[X] — Et^x ^pc = Et^x 



^PC=Et,x PodX{T) ~h{X{T),T)dT (4.10) 



where h is now a scalar field, called the energy function. 
For (j4.ip it reduces to 

MX, r) = + • (D.X) ViX) (4.11) 

We shall come back later on the meaning of h. 

Assume now, following Feynman, that we are only given the Lagrangian C 
of (|4.7p . For the sake of generality we shall add a smooth (final) boundary 
condition: J[X] = Ex,t{J^CdT + Su{X{u))}. 

When not referring to the C of (|4.ip we shall assume from now on that 
C,Su are continuous and that , for some constants c,k, \£\ ^ c(l + \X\'') + 
\DX\MS^iX)\^cil + \X\'^). 

How can we characterize critical points of J ? 

We shall pick, as domain Vj of J the set of diffusions X absolutely contin- 
uous with respect to the Wiener measure with diffusion matrix hi (J the 
3x3 Identity matrix), and arbitrary Borel measurable drift B oiM.^ x [s, u] into 
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Def. Such a process X is extremal for J if 



E,t[VJ[X]{SX)]=E,t 



',. J\X + eSX] - J\X] 
lim 



0, 



(4.12) 



for any variation 5X in the Cameron- Martin space preserving the absolute con- 
tinuity under the shift X + eSX (cf fl^, fE^, |2f ). 

Then 



dr- 



E,t[VSu{X{u))dX{u)]. 

In the second term, notice that we preserve the notation D^, this time for the 
(extended) infinitesimal generator of X G Vj, 



d h , 

i?, = — + SV + -A. 
OT 2 



(4.13) 



When apphed, in particular, to X itself, Dt-X — B{X,t), as in (I3.15p . Under 
Ext[- ■ ■], D-r satisfies, by Ito formula, an integration by parts formula. Since 
5X[t) is of bounded variation, 



dC 
dX 



Dr 



dC 
dDrX 



SX{T)dT+ 

f dC 



E. 



\dDrX 



+ V5„ {X{u))5X{u) 



so we obtain the 



Stochastic Euler-Lagrange Theorem 

The critical processes of the Action functional J[X] with boundary condition 
Su solve the almost sure Stochastic Euler-Lagrange equation: 

[Dr[s^)-§^^Q. t<r<u 
(SEL) <^ (4.14) 

i DrX{u)) = -VSuiX{u)), X{t) - X. 

For instance, when C is as in (14. ip . (SEL) reduces to 



DrDrX{T) = VV{X{t)) + DrX{T) A rot A + f rot rot A 
DrX{u) + A{X{u)) = -ySu[X[u)). 



(4.15) 



One should notice the h dependent term, in addition to the Lorentz force, on 
the r.h.s of Eq (|4.15p . It is natural to interpret (SEL) as the {Vt) stochastic 
deformation of its classical counterpart. 
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Regarding the Hamiltonian approach of the dynamics, one can define, as 
classicahy, using the hypothesis after (|4.8p . 



(4.17) 



H{X, P) = P</.(P, X) - C{X, </.(P, X)). (4.16) 

The deformed version of the Hamiltonian differential equations becomes the 
almost sure 

Stochastic Hamiltonian equations: 

iDrX^^ 

For instance when C is the one of (|4.ip . H is 

mx, P) = l\P- AiX)\' - V{X) - • A (4.18) 
and the Hamiltonian equations reduce to 

(SHE) 




DrX A rot A + frot rot A+ ^AA + W. 

(4.19) 

The first equation reduces to another version of p.lSp . Its substitution in the 
second one is consistent with the Euler-Lagrange equation (|4.15l) . 

Almost sure equations like SEL (|4.14p and SHE (|4.17p may seem to be odd 
but, in fact, they were (rather deeply) hidden behind classical results of Stochas- 
tic Control Theory found around the eighties [T8j [T9j [2^ from a very different 
viewpoint. In this context, it is the classical theory of Characteristics which was 
stochastically deformed. We shall give here only a hint about this connection, 
particularly natural in our Path Integral perspective, for the simplest case above 
where A = 0. 

Consider, then, the Action functional J[X], where the drift B of the critical 
process (with fixed diffusion matrix) we are looking for, called the "control" is 
just supposed to be Pt-measurable and such that E |i?(r)|"(iT < oo, n € N. 
Notice that this include, now, non Markovian processes, a natural hypothesis in 
our Bernsteinian perspective. 

Consider any scalar field S in the domain of the infinitesimal generator 
A^^"^) = + S(r)V + |A of such a process X{t), such that 

Ext J \ A^^'''> S {Xr, t)\ dr <oo 

and Dynkin formula holds: 

Ext ^ A^(^)5(X,, r) dT = ExtS{X{u), u) - Six, t). 
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So, for any Xlj) in the above class, and £(X, DrX) = \\DrX\'^ + V{X) 

J[X] = E^t ^ C{X{t),B{t)) dT + E,tSu{X{u)) 

Theorem ([18'). Let S{x,t)be a classical solution of the deformed Hamilton- 
Jacobi equation (known as "Hamilton- Jacobi Bellmann", or HJB) 

r|2_i|V5P + |A5 + y^0 ^^^^^ 

ys{x,u) = Su{x). 

Then S{x,t) ^ ^[-'^l, VX in the above class. Moreover, for 

B{t) = B{x, t) = -yS{x, t) (4.21) 
this inequality becomes an equality. 

To understand the relation with our construction (for A = 0) , define 

77(3;, t) ^e-*'^^"'*) (4.22) 
then HJB reduces to Eq (|3.1ip with a positive final condition, 

ft^ = -yA77 + y, 77(x,?/) =e-*^"(") (4.23) 

and the critical, indeed minimal, diffusion of J[X] is our Markovian (forward) 
Bernstein process of drift p.lSp DtX = B{X,t) = h\7 log ri{x,t). 

It is in fact possible to prove in a purely geometric way that the critical 
points of J[X] are minimal pT| . 

The hypothesis, in the last Theorem, that 5 is a classical solution of HJB is, 
of course, much too restrictive. The construction holds under weaker conditions 
(cf. [H]). Regarding HJB, the appropriate notion of weak solution is the one 
of viscosity solution (cf. [IB]). 

When h = 0, the existence of a classical, global solution of Hamilton- Jacobi 
equation is a condition of Complete Integrability of the system and, if available, 
the gradient of this equation coincides with the second Hamiltonian equation, 
f = -It = -Vl^(x), when nix,p) = + vix). 

The stochastic deformation of this integrability condition is that, the gradi- 
ent of (I4.20p . using DtX = W log 774 (i.e p.lSp for A = 0) reduces to, almost 
surely, DtDtX{r) — W{X{t)) namely (|4.15p or, equivalently, the second de- 
formed Hamiltonian equation (|4.17p in the same special case (cf. [23]). 

Coming back to the deformed Hamiltonian H of (|4.18p ones observes that 
it does not coincides with the energy function h (|4.1ip of the Poincare-Cartan 
1-form. So, what is the meaning, if any, of h{X{T),T) 1 

A key observation is that Drh{X{T),T) = 0, for the critical, i.e dynamical, 
diffusion X{t), with generator 

= A + ftviog77V + Ja. (4.24) 
or 2 
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In other words, /i(X(t),t) is a 'Px-martingale, a natural deformation of the 
classical notion of constant of the motion or first integral. Is it accidental ? 

The answer is negative. Before mentioning the Theorem showing that all 
first integrals of our stochastic dynamical system are indeed martingales, let us 
stress that Feynman could not find such a result in his informal time discretized 
account of quantum dynamics. 

As it is clear from (I3.15[) or (|4.19p . the processes critical for J are entirely 
built from (positive) solutions of a parabolic equation p. lip . We shall stick to 
the simple case A = 0, i.e Eq (|4.23p . 

Consider this equation, written now as 

f d \ 
Hri^ [h— + —A - 1/ 7/ = (4.25) 



at 2 



and a generator of the form 



o (5 1 

N^Tit)- + Q,ix,t)—--cbix,t) (4.26) 

where the summation convention is used, for 1 ^ j ^ 3 and the coefhcient 
T,Q,(j) are analytic in x and t. Then N is infinitesimal generator of a symmetry 
Lie Algebra As of Eq (IT^ if 

Hr] = 0^ HNt] = 0. (4.27) 

Lie proved long ago that this is the case if (T, Q, (f) solve the following "deter- 
mining equations" : 

'^=2^ |2. + ^=0, 7-1,2,3, 

at oxi oxj oxi ' -! 1 1 J / 

dQj_ ^ dip 
dt dx^ 

If + |A0 = ^T^ + Q.C + r^ 

(4.28) 

The (local) Lie symmetry group of Eq (|4.23p results from product of expo- 
nentials of generators N. 

When a; € K'^ and V — A — 0, ioi instance, the symmetry group is 13 
dimensional. 

Now consider a classical Lagrangian L{uj, uj, t) associated with the deformed 
system whose Eq (|3.1ip reduces to the form (|4.23p . it is L{uj,di,t) = ^I^P — 
V{uj,t) (The scalar potential V can be explicitly time-dependent without alter- 
ing our arguments above) . One version of the classical Theorem of Noether is 
the following: 

d ^ . . d fdQi . dT\ d 



T{t)-+Q,{uj,t) 



dt ' duji \ dt dt J duji 

is called a divergence symmetry of L if there is a scalar field (p = '/'(w,t) such 
that 

/ , N ^ dT d6 , , 
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When L admits such a divergence symmetry, then along each extremal of the 
classical action Sl, 



d 
H 
d_ 
di 



—I 

duji 



dL 

dio. 



■u},, - L ] T - 



(4.30) 



'^^Q^- ( ^\u;\^ + Viuj)]T^cj) 



= 



for L as before. Equivalently, the expression between brackets is a constant of 
the motion of the system. 

For instance, any conservative system (i.e with V time independent) admits 
V = -g^ i.e T — I, Q ^ 0, (p ^ and (|4.30p reduces to the energy conservation. 

The stochastic deformation of this Theorem, for the same class of systems 
is the 



Stochastic Noether's Theorem [25l[26] 

If the Lagrangian C{X, DtX, t) = ^\DtX\^ + V{X^ t) admits a divergence sym- 
metry of the form 

for any analytic T,Q,(p solving the "Determining equations" (I4.28P then along 
any Bernstein diffusion X{-) critical for the action J[X], almost surely 

Dt{DtX,Q,~hT-q^){X{t),t)=0 (4.32) 

where h is the energy function (|4.1ip . 

For instance, if V is time independent, C admits T — 1, Q^(l)^0 and we 
recover our energy martingale. 

The stochastic Noether Theorem is a theorem of structure, here, without 
which our deformation would be dynamically meaningless. 

But let us observe that, from the start of this Section 4 only one filtration, the 
increasing one Vt, s < t < u, has been used. As a result of this, the stochastic 
Euler-Lagrange equation (|4.15p . for instance, is certainly not invariant under 
time reversal in the sense defined in Section 3. So (|4.15p cannot be the full 
dynamical characterization of processes respecting intrinsically this invariance. 

But it is easy to find the solution of the puzzle. On the time interval [s, u], 
the time reversal of the Action functional (|4.7p is 

E^' J' (^\\D;Xir)\' + V{Xir))^ dr + j\odX{T) (4.7*) 

where we have adopted the notation E^^ for a conditional expectation given the 
future configuration X{t) — x, s <t, and used the rule p.23p . 

Calling J* [X] the functional ()4.7*P , we can look for its critical points among 
diffusions with fixed diffusion matrix, solving a backward (i.e Tt) stochastic 
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differential equation whose (backward) drift is unknown. Notice that, now, the 
Stratonovich integral in ()4.7*l) must be interpreted using the backward version 
of the relation (|4.6p namely, according to Ito [13] , 

AodX{T) ^ Ad^X{T) - ■ AdT. (4.6*) 
This means that, with respect to J>, the Lagrangian of J* is now represented 

by 

C*iXir),D:X{r)) = l\D*X{r)\^ + V{XiT)) + 

AiXir))D;X{r)^^V-A. (4.1*) 

Then, one checks in the same way as before, that the critical point of J*[X], 
in fact a minimum, is unique and that its backward drift is given in term of a 
positive solution of Eq (|3.10p by the expression ()3.15*|) : 

D;X = -h\7\og7j;{X) - A{X). 

The (backward) stochastic Euler-Lagrange equation it solves (ignoring boundary 
condition at i > s), is 

D*D;X{t) = (X(r)) + D*X{t) a rot a - ^rot rot A. (4.15*) 

As a matter of fact, such a calculation is not even necessary. Indeed, as said be- 
fore (cf (|3.23p '). Dr — s> —D* and A — > —A under time reversal. This means that 
both (j3.15*|) and (|4.15*|) are time reversed versions of their forward counterparts 
((31^ and (|4J5l) . 

Since, in particular, ()4.15p and ()4.15*|) provide different informations, as- 
sociated with Vt and J-t respectively, about the same Bernstein diffusion, the 
complete, time-symmetric, dynamical equation of X{t)^ s ^ r ^ m, is 

^{D,DrX{T)+D;D;X{r)) = 

^{DrX{T) + D*X{t)) a loi A{X{t)) + VV{X{t)). (4.33) 

Let us stress that, now, this stochastic deformation of the classical Euler- 
Lagrange equation is an electromagnetic field, involving the deformed Lorentz 
force on the right hand side, is indeed invariant under time reversal, as it should. 
Using the relations , for / G C^, 

l-E[f{X{T))] ^ E[D,f{X{T))] = E[D*J{X{t))] 

following, for instance, from Dynkin formula (cf. also ^53j), we can get closer to 
Feynman dynamical law (|2.18p in taking the absolute expectation of (|4.33p : 



^.ElXir)]^E 



\{D,X{t) + D*X{t)) a rot A(X(r)) + VV{X{t)) 



(4.34) 
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Notice also that what plays the role of the time derivative uj{t) in Feynman's 
law of motion (|2.18l) is now the average of the two drifts, whose behaviour 
under time reversal is the expected one, in contrast with each of them taken 
separately. Using the same method, it is easy to find the backward version of 
our Stochastic Noether Theorem, for instance, producing backward martingales 
of the system. 

5 Computational and Geometric content 

Let us start with some consequences of our Noether Theorem. 

Although it is clear, for a member of the community of Geometric Mechanics 
(in particular) that Noether is the key to start a serious study of the dynamics, 
we shall try to show its interest also for the theory of stochastic processes itself. 

Consider diffusions on the line, for simplicity, with A = V = in Eq p.l3p 
or, equivalently a Lagrangian (|4.1I) reduced to 

C{X,DtX) = ^\DtX\^. (5.1) 

One verifies that T — 2t, Q ~ x, (j) ~ solves the one dimensional version of 
the determining equations (14.281) for = 0. So N = 2t-^ + x-^ generates a 
one-parameter symmetry group: 

(e"^)(i,x,77) = (e"x,e2"t,ry) = (i,, ?7a). (5.2) 

This implies that if r] = ri{x, t) > solves the free heat equation (|4.23p for 
y = 0, i.e h^ = -'4^,8o does 

rj^{x, t) = e-"^77 = 77(e-"="x, e"^"i)- (5-3) 
Then define K{x,t) = '-^{x,t). If X{t) solves the (forward) SDE with drift 

dX(t) = nV\ogr]{X{t),t)dt + h^^^dWt (5.4) 
and, therefore, the associated a.s Euler-Lagrange Eq (I4.15|) . 

DtDtX{t) = (5.5) 

one checks that ha{X{t),t) is a positive martingale: 

Dtha{X{t),t)^0. (5.6) 

Eq (|5.3p corresponds to a Scaling transformation of the starting process X{t), 
namely 

X"(t) = Xo.{X{t^),t{ta,)) = e"X(e-2"t) (5.7) 
where t{ta) denotes the inversion of the time parameter transformation in (j5.2p . 
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The drift of X^it) results of a Doob's h transform of X(t) whose martingale 
is ha- Denoting by B" and B the associated drifts, we find [27], using the 
definition of ha, 



In other words, all diffusions X{t) solving Eq (j5.5|) enjoy the scaling transforma- 
tion symmetry, interpreted here dynamically. The standard Wiener X{t) = Wt, 
whose ria{x,t) = 1, therefore ha = I and Ba = 0, is only one of them: with 
e = e~^", we recover its usual scaling law W^^(t) = e~^^'^W{et). This is useful, 
for instance, in the computation of first passage times of any diffusions solving 
Eq (|5.5|) . A large collection of parabolic equations (with first order and poten- 
tial terms) is, in fact, equivalent to the above free heat equation, in terms of 
symmetries [27], so that our argument is more general. 

Almost sure dynamical equations like (|4.14l) or (j4.17p . together with our 
stochastic Noether Theorem provide new (geometrical) relations between famil- 
iar stochastic processes, impossible to anticipate without them. But, what may 
seem more surprising, they provide as well new informations about Quantum 
Mechanics in Hilbert space. For instance, the solutions (T, Q, (p) = (0, t, x) of 
the (one dimensional) determining equations (|4.28p for y = 0, applied to the 
standard Wiener X{t) = Wt, show that the family of Brownian martingales 
correspond to a family of (time dependent) constant observables in L^(R) of 
the free particle. Cf. [3^. Although elementary, this observation had not been 
done before. Even better, a naive analytic continuation in time from the sym- 
metries of the parabolic equation to the one of Schrodinger equation provides 
a Quantum Theorem of Noether richer that the one mentioned in Textbooks, 
even in elementary cases [29j . 

The geometrical content of our stochastic deformation is worth an investi- 
gation in itself. Consider, for instance, the deformation of the classical method 
of Characteristics. One of the most elegant representation of Hamilton- Jacobi 
(HJ) equation is due to E. Cartan and makes use of the following Ideal of dif- 
ferential forms [3^ (We denote here by Xi what was w^, in Eq (|4.9p . to avoid 
confusions with differential forms): 



on a 9-dimensional space {xi,p^ , S,t, E) i = 1,2,3. 

This is Cartan's representation for the elementary systems underlying Eq 
(|4.30p . To recover HJ equation itself, consider the E'* "solution submanifold" 
where the a priori independent variable S becomes a function S{x, t) (This is 
called "Sectioning" and denoted by ^) and then puUback all differential forms 
to zero ("Annuhng"). Then w = implies that p = -VS* and E = drS. The 
condition = is equivalent to the existence of a Lagrangian manifold. Finally, 
/3 = is equivalent to the classical Hamilton- Jacobi itself. This representation 



B°'{x,t) = B{x,t) + V log ha{x,t) = h\/ log r]a{x,t). 



(5.8) 




(5.9) 
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shows that Hamilton- Jacobi framework is a Contact Geometry, defined on an 
odd dimensional space, here MP; ui is, in fact, a Contact form |27j . 

A symmetry generator of HJ equation should, therefore, became now a "con- 
tact Hamiltonian" vector field (sometimes called "Isovector") 

whose coefficients must be chosen so that, denoting by Cn the Lie derivative, 
or variation, along N: 

Cn{Ih.j) C Ihj. (5.11) 

The stochastic deformation of Ihj is the one providing Hamilton- Jacobi Bellman 
equation (14. 20^ : 

{uj = P^dX, + EdT + dS = ujpc + dS 
n = dP'dX,+dEdT (5.12) 
I3 = {E + i|P|2 - V) dXidr + ^dP'dr. 

The only deformation term, in /3, is responsible for the deformation term fAS* 
in Eq (I4.20p . Sectioning and annulling as before we find the Lagrangian inte- 
grability conditions: 

= 0^ P = -VS,E = ~drS. (5.13) 

The definition of Symmetries for HJB is the same as classically, i.e Eq (|5.11[) . 
for the deformed ideal (|5.12p . and the calculation of the coefficients N* of Eq 
(|5.10p is a tiring exercise (cf. |27|). But it is quite rewarding: 

Theorem ([27]). Along any N -variation as before, Irjb o,nd the Lagrangian 
C satisfy the following invariance conditions 

(1) CNii^pc) = -dN^ 

(2) Cn{^) = (5.14) 

(3) Cn{C) + = -DrN^ . 

dr 

This Theorem seems purely algebraic but encodes a lot of informations about 
our stochastic deformation, resulting from the substitution of smooth classical 
paths T i-T- uj{t) by Bernstein diffusion sample paths r i— > X(t). Eq (1) means 
that Poincare-Cartan 1-form is invariant up to a phase coefficient N^. Eq 
(2) shows the invariance of the Symplectic form over the time-dependent or 
extended phase space (cotangent bundle). Eq (3) expresses the transformation 
of the integrant of the Action functional (|4.7p under the contact Hamiltonian A'' 
on the extended phase space. It should be regarded as the deformation of the 
classical expression (|4.29p . 

The proof of the Theorem shows that it is, in fact, sufficient to consider 
symmetry contact Hamiltonian of the form N{t,x, S, E, P) — N''{x,t)P + 
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N^{x,t)E + N^{x,t), so that N"" ^Q, = T and = in the notations 
of the stochastic Noether Theorem, where T, Q and (p solve its Determining 
Equations. After sectioning on the solution submanifold (x, t) where, according 
to p.9|) . we have a probabilistic interpretation by plugging x = X{t), P and E 
become respectively our drift and energy random variables. 

The whole construction summarized before is preserved if the diffusions X{t) 
lives on a (smooth, connected, complete) n-dimensional Riemannian manifold 
with metric gij. The simplest Hamiltonian in Eq (j4.25p becomes 

^ = -yV^V, + y(x) (5.15) 

where Vj denotes the covariant derivative with respect to the Levi-Civita con- 
nection r*^, and the associated Hamilton- Jacobi-Bellman equation (|4.20p be- 
comes 

^~l\\^Sf + ^yV,.S + V^O. (5.16) 

Two new geometric aspects deserve to be mentioned. The first one is that (as 
stressed by K. Ito [31 ) an additional term shows up in the drift p.l5p 

DtX^ = h\7 \og7n{X) - ^V)^g^^ (5.17) 

for EtdX'dX^ = hg'idt. 

On such a Riemannian manifold, an almost sure Euler-Lagrange equation 
like ()4.15p (when A — Q) requires to define the time derivative of a vector field. 
Even in the classical, deterministic case, a notion of parallel transport is needed 
to do that. According to Ito [JT], the stochastic deformation of the Levi-Civita 
transport of the vector field Y would transform Eq (|4.13p into 

DrY' = — + hVHogiitVkY' + -V'=Vfcr\ (5.18) 
OT 2 

But Ito also indicated other possible choices. The one needed for our purpose 
has been called "Damped parallel transport" in Stochastic Analysis (cf. ^35]), 
and replaces the Laplace-Beltrami term of (I5.18P by 

(Ar)* = V'^VfeF^ -f WkY'' (5.19) 

where R^. denotes the Ricci tensor. Then 

fiyi ft 

DrY' - ^ + W^logr^tVfey' + -(Ay)\ (5.20) 
OT 2 

The point is that to preserve for ()5.16p the integrability condition according 
to which the gradient of (|4.20p coincides with the Euler-Lagrange equation we 
need that [A, \I^]S = 0, a property not satisfied by V'^Vfc. Then, with (|5.20p . 
the dynamical law and the Noether Theorem keep the same form as above |61| . 
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6 Conclusions 



In 1985-6, I named after Bernstein the reciprocal property suggested by him 
in the context of the 1931 observation (forgotten until then) of Schrodinger 
[B]. I was, in fact, so impressed by Bernstein interpretation of such processes as 
stochastic counterparts of critical trajectories of Hamilton's principle that I used 
as well the term "variational processes" P]. Of course, the local Markov prop- 
erty reappeared during the seventies in relation with Statistical and Quantum 
Physics. But Schrodinger's observation and Bernstein's probabilistic sugges- 
tion were an extraordinary anticipation back then, of Feynman Path Integral 
approach, among other ideas. 

There are more than one way to interpret Schrodinger original observation, 
expressed originally in a statistical mechanics perspective, i.e in entropic terms. 
Aware of the fierce fights regarding the physical interpretation of the new born 
quantum theory, Schrodinger was looking for a "classical" analogy where prob- 
abilities would play a similar but less debated role. In the seventies, B. Jamison 
[8] elaborated some aspects of the construction suggested by Bernstein, but 
missed the looked for relation with quantum theory (he was, for instance, using 
only the increasing Vt filtration and would not start his construction from a 
given Hamiltonian H). In any case, since the mid-eighties, Bernstein processes 
have reappeared in a multitude of contexts, pure and applied, and under differ- 
ent names. They have been called "Schrodinger processes" (following Jamison) 
by H. Follmer |33j in 1988 and studied on their own in the entropic perspective 
[33], [in] and [50] ■ A promising link has been established with Optimal 

transport in recent years. An excellent review of this connection can be found 
in |36j . In this context, the natural approach is indeed the one of statistical 
mechanics and the original variational problem in [B] is called "Schrodinger 's 
problem", not to be confused with Eq p.lSp referred to as the "Schrodinger 
system". "Schrodinger bridges" is also a terminology used for those processes. 
They can really be regarded as a generalization of usual bridges where, instead 
of two boundary Dirac distributions, we are now given two regular (nodeless) 
probability distributions. In recent studies on Wiener space, they have also 
proved to be quite natural tools [37] . 

Reciprocal Bernstein processes can also be characterized by an integration 
by parts formula, typical of Stochastic Analysis, but even when they are not 
Markovian [38) . It seems, indeed, that Feynman's symbolic approach was too 
limited to the Markovian class, not appropriate in some cases. 

As it is clear from the first and last part of Section 4, we need 3 kinds of 
stochastic integrals for the complete description of Bernstein processes. A very 
general approach to stochastic integration with the same features (besides the 
one of ^.) due to Russo and Vallois is known: of [60]. It would be interesting 
to reconsider our construction with the tools described there. 

The symmetry of such processes, in the sense of Noether Theorem, can also 
be of interest for other purposes [Slj • 

Clearly, the approach chosen here can be regarded as a random version of 
Geometric Mechanics (cf. [23], [54]). In this context, one of the most interesting 



26 



open problem is the notion of Integrability suited to the random dynamical 
systems resulting from our approach. Some aspects of it have been used in the 
Ideal of differential forms (j5.12p but a lot more remains to be done (cf. [S^) in 
this field. 

It would be interesting as well to understand the relations between our 
stochastic deformation and the (deterministic) deformation of characteristics 
for Hamiltonian PDEs inspired by B. Dubrovin [62]. One can, indeed, guess the 
existence of common features, for some particular PDEs. 

Some probabilists would wonder why to study such a special class of Bern- 
stein reciprocal processes. The first reason is that this class is not as small as it 
seems. We hope that we made clear that the key elements of their construction 
are independent of the form of the starting Hamiltonian H . Besides those like 
(|3.25l) . for Levy processes, we claim that it is always possible to time-symmetrize 
regular stochastic processes the way we did here. 

The second reason is that if it is precisely because this class is special that it 
carries all the qualitative properties needed to construct stochastic dynamical 
theories. As mentioned in §2, Stochastic Analysis did not go at all in this 
direction. But, as suggested by Feynman Path Integral approach, this direction 
seems to be the most natural one as far as physical theories are concerned. 

There are many fields, outside Mathematics, where this unorthodox way to 
approach stochastic dynamics is also natural. For instance, in Finance [32] or 
Econometry |40| . Image processing is also a promising domain of application of 
these ideas. 

Of course, various problems of Statistical Physics can benefit from the use 
of such time reversible probability measures [41j . P.O. Kazinski, among many 
others, considers various classical models in this perspective [42]. He is also 
the author who introduced the expression "Stochastic Deformation". Other 
applications in Theoretical Physics include [43], [58] . 

Random walks on graphs are described by Markov chains, reversible in a 
narrower sense than the one intrinsic to Bernstein processes. It is likely that 
the methods used here will also be relevant in this area [S^. Interesting links 
with physics are explored in [5^ . 

In Applied Mathematics, the relations of the variational component of the 
program summarized here and Stochastic Control Theory are, of course, strik- 
ing. They strongly suggest that there are very few ways to deform systematically 
classical mechanics along diffusion processes. But those relations are still far 
from completely explored. It is remarkable, as mentioned before, that some in- 
vestigations of the seventies and eighties, aiming at a deformation of the classical 
calculus of variations along diffusion processes, were able to obtain results con- 
sistent with our probabilistic reinterpretation of Feynman's approach. What 
U.G. Haussmann |20j calls the adjoint process, for instance, is basically our 
(forward) momentum process (|4.8|) . Of course, those results were all expressed 
with respect to a single (increasing) filtration and, as such, were not directly 
appropriate to a time reversible dynamical framework. 

A last comment about the Stochastic Deformation program. In the late 
sixties, V. Arnold proved that the Euler equation of an ideal incompressible 
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fluid could be interpreted in a (Lagrangian and Hamiltonian) analogy with the 
motion of a rigid body [44]. The configuration space was, then, the group of 
volume preserving diffeomorphisms of the region occupied by the fluid. If we are 
not only interested by "dry water" (as Feynman was calling the fluid described 
by Euler equation [55]) then we have to deal with Navier-Stokes equation. The 
idea that Navier-Stokes equation corresponds to a stochastic deformation of the 
Euler one was introduced in [52) and has been considerably elaborated in recent 
years [5S] . 

This means that the method of stochastic deformation can also be applied 
to some infinite dimensional dynamical systems. 

Various books have already been published, where Bernstein reciprocal pro- 
cesses play a major role. We mention only two recent ones besides [28]: [47] . 

m- 

We are convinced that, also on the applied side, those processes do not have 
only a curious past but also, indeed, a bright future. 
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